Detecting nonclassical system-environment correlations by local operations 
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We develop a general strategy for the detection of nonclassical system-environment correlations 
in the initial states of an open quantum system. The method employs a dephasing map which 
operates locally on the open system and leads to an experimentally accessible witness for genuine 
quantum correlations, measuring the Hilbert-Schmidt distance between pairs of open system states. 
We further derive the expectation value of the witness for various random matrix ensembles modeling 
generic features of complex quantum systems. This expectation value is shown to be proportional 
to a measure for the quantum discord which reduces to the concurrence for pure initial states. 
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In the theoretical analysis of the dynamics of open 
quantum systems it is often assumed that the open sys- 
tem S and its environment E are statistically indepen- 
dent at the initial time. Under this assumption the re- 
duced dynamics of the open system can be described in 
terms of a completely positive dynamical map which op- 
erates on the open system's state space, mapping the 
states of S at time t = to the states of S at time t > 0. 
The one-parameter family of such maps indexed by the 
time-parameter t then allows the description of Marko- 
vian as well as non-Markovian quantum processes (H-Q. 

The mathematical formalism of completely positive 
transformations provides a fundamental concept in the 
quantum theory of open systems [3( and in quantum in- 
formation theory [Bj]. However, the assumption of an ini- 
tially uncorrelated state is in general physically hard to 
motivate. The description of the reduced system dynam- 
ics in terms of linear completely positive maps fails in the 
case of general initial correlations. Several approaches 
have been proposed to treat the problem of the dynam- 
ics for initially correlated states (see, e.g., j(i-[l6j]). 

One of the central problems in the treatment of cor- 
related systems is given by the fact that the total initial 
state is typically unknown since a full experimental con- 
trol of the environment is either impossible or much too 
complicated in practice, while only the reduced subsys- 
tem is experimentally accessible through local measure- 
ments and operations. The present work is thus moti- 
vated by the question, how can one detect the presence 
of initial correlations in an unknown system-environment 
state p through only local operations carried out on the 
open system? Our approach builds on the strategy pro- 
posed recently in Ref. [17] . where it is shown that cor- 
relations in two total initial states p and p' can be wit- 
nessed by comparison of the time evolution of their re- 
duced open system states ps(t) and p'gif), respectively. 
Very recently, experimental realizations of this scheme 
have been reported [3, [l9[ . In order to witness corre- 
lations in a given unknown state p, the basic idea is to 
create the second state p 1 — ($ <g> I)p by applying a local 
trace preserving operation $ to the original state p. Since 



a local operation can never generate system-environment 
correlations, any increase of the trace distance between 
the open system states over its initial value implies the 
presence of correlations in p [171 ] . 

The above scheme may fail to detect initial correla- 
tions. Whether or not it is successful depends crucially 
on the chosen local operation and on the structure of the 
unitary time evolution operator Ut which describes the 
total system dynamics. Moreover, the general scheme 
does not explicitly distinguish between classical correla- 
tions and genuine quantum correlations of system and 
environment. The central goals of the present paper are 
thus twofold. First, we develop a general strategy which 
allows the detection of true nonclassical correlations, i.e., 
of correlations described by states with nonzero quantum 
discord. This is achieved by taking the local quantum op- 
eration $ to be a dephasing map describing the complete 
decoherence in a suitable local basis. The amount to 
which the evolution of the reduced system states ps(t) 
and p'g(t) differ from each other can be quantified e.g. 
by means of the Hilbert-Schmidt distance between these 
states. Our second goal is to investigate the generic be- 
havior of this distance, without making specific assump- 
tions about the system and its total Hamiltonian. To this 
end, we determine the expectation value of the Hilbert- 
Schmidt distance for various random matrix ensembles 
describing typical features of complex quantum systems. 

We consider a state p of the total system and its re- 
duced open system state ps = Txe/) given by the partial 
trace over the environment. While the total system state 
p is typically unknown, we assume that the reduced state 
ps can be determined through state tomography yielding 
a complete set of rank-one projection operators 7r^ which 
project onto the eigenvectors of ps- Employing these op- 
erators we define a local dephasing operation $ acting on 
operators A of the open system by means of 



(1) 



Note that the 7r M and the quantum operation $ depend 
on the total system state p, but the full information nec- 
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essary to construct these quantities is locally accessible 
through measurements on the open system. The appli- 
cation of the map $ corresponds to a non-selective mea- 
surement with respect to the basis given by the projectors 
7r M , and can also be interpreted as describing the com- 
plete decoherence of ps in this basis. The corresponding 
change of the total system state is given by 



P 



' = ($®/) P = ^ivn A 



(2) 



where / denotes the identity map and II M = tt^ (g) /. 

By construction, the marginal states of system and en- 
vironment are invariant under the application of ($(8)/), 
i.e., we have p' s = Tr^p' = ps and p' E = Tisp' = Pe- 
However, in general the total system state does change 
under the local dephasing operation. The requirement 
p = p' for the invariance of the total state under the 
dephasing operation can be written as 



(3) 



This equation represents the condition for the total sys- 
tem state p to have zero quantum discord, i.e., to contain 
only classical correlations [2(| 21 1 . The quantum discord 



has been introduced as an information theoretic measure 
for the quantumness of correlations in a composite state 
p. The original definition involves two different quan- 
tum generalizations of the classical mutual information. 
Here, we use a technically simpler measure for the degree 
of quantum correlations as suggested in 22j : 



6[p) = \\p-{*®I)p\\, 



(4) 



where \\A\\ = VTrA^A denotes the Hilbert-Schmidt norm 
of an operator A. The quantity Q provides a measure for 
the change of the given state p induced by the measure- 
ment of the local eigenbasis of the reduced open system 
state. It vanishes if and only if Eq. ([3]) holds true, which 
means that there exists a local measurement basis which 
does not disturb the total system state, motivating the 
notion of classical correlations. It is useful to note that 
the expression Q can also be written in terms of the 
purities of the states p and p' . Defining the purity as 
V(p) = Tr{p 2 } we have 



6 2 (p)=V(p)-P(($®I)p), 



(5) 



which can easily be verified by use of the relation Tr{ [($(£> 
I)p] 2 } = Tr{p($ ® I)p} which in turn follows from ((2|). 
Thus we see that the measure (j4|) for quantum corre- 
lations is the greater the more the purity is decreased 
through the local dephasing operation. 

Consider now the dynamics of the total system given 
by some unitary operator U t on the total system's Hilbert 
space H. We denote by ps(t) — Tr~E{U t pU}} and 



p' s (t) — TiE{Utp'Ut} the open system states correspond- 
ing to the total initial states p and p' = ($ ® I)p, respec- 
tively. Although the open system states are identical at 
the initial time t = 0, the time-evolved states ps(t) and 
p' s {t) can differ at some later time t > due to different 
correlations in the total initial states. This means that 
different correlations in the initial states become dynam- 
ically perceivable in the open system states. If ps{t) and 
p's(t) evolve in fact differently one can conclude that the 
initial state p must contain nonclassical correlations since 
otherwise p and p' are identical. Thus, we find that the 
condition 



Tr E {u t (p-p')U}} 



> 



(6) 



implies the presence of nonclassical correlations in the 
total state p, i.e., S(p) > 0. 

To illustrate the application of our scheme we refer to 
a recent experiment [19( in which initial correlations be- 
tween the polarization degree of freedom of a photon, rep- 
resenting the open system, and its translational (mode) 
degree of freedom, forming the environment, have been 
observed. To detect initial quantum correlations in this 
experiment one first performs a polarization state tomog- 
raphy on the photons leaving the optical fiber which is 
used in the experiment to create a correlated initial state. 
This allows to determine the eigenbasis of the photonic 
polarization state. The action of the local dephasing map 
4> can then be realized by placing appropriate polariza- 
tion filters right after the optical fiber (before the pho- 
tons enter the delay setup), in order to carry out a non- 
selective measurement of the polarization state in this 
eigenbasis. The quantity ^ can then be measured by 
state tomography after the delay setup. 

The condition yields a general method for the de- 
tection of nonclassical correlations by local quantum op- 
erations. However, this condition is only sufficient for the 
existence of nonclassical correlations since it may hap- 
pen that a particular Ut does not lead to an increase of 
the Hilbert-Schmidt distance in the reduced state space, 
even though p is nonclassically correlated. The question 
is thus, can one derive general statements about the be- 
havior and the size of the distance in ^ for a generic 
choice of the time evolution operator? To answer this 
question we regard the time evolution operator as a uni- 
tary matrix U which has been drawn from an appropri- 
ate random matrix ensemble over the unitary group U(d), 
where d = dsds with d$ and rig denoting the dimensions 
of the Hilbert spaces of the open system S and its envi- 
ronment E, respectively. The expectation value of some 
functional F(U) will be denoted by angular brackets, 



(F(U)) ee / dp(U)F(U), 



(7) 



where dp(U) represents the probability measure on the 
unitary group U(d). As a first step we consider the ma- 
trix ensemble given by the uniform Haar measure oiU(d). 
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Our goal is to determine the average Hilbert-Schmidt dis- 
tance between the open system states with respect to the 
Haar measure. 

Theorem. Let M be an arbitrary fixed Hermitian 
operator on H and define the reduced system operator 
A = Tr E {UMW}. Then we have: 
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(8) 



The proof is given at the end of the paper. To apply the 
theorem to the present case we set M = p — p 1 which is 
an operator with trace zero. Therefore, the second term 
of Eq. JHJ vanishes and we find: 



\T?v{U{ P - P ')W}\\ 
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Sip). (9) 



Thus we obtain the remarkable result that the root mean 
square of the Hilbert-Schmidt distance between the open 
system states is proportional to the measure (j4]) for non- 
classical correlations in the initial state p, with the prefac- 
tor depending only on the dimensions dg and d E - While 
the witness ([6]) may fail to detect initial quantum correla- 
tions for a particular unitary evolution operator, Eq. (j9]) 
shows that the expectation value of this witness with re- 
spect to randomly drawn unitaries is nonzero if and only 
if the initial state contains nonclassical correlations. 

An important special case of the general result (HJ is 
obtained if p = | 1 J r )(^ , | is a pure state, such that entangle- 
ment is the only possible type of correlations. Employ- 
ing the Schmidt-decomposition 1^) = J2 i X i \ipi) £g) 
it is easy to see that the local dephasing operation 
projects onto the local Schmidt basis |</jj), which yields 
p' = Y,t A f \<Pi)(<Pi\ ® \Xi)(Xi\- Using again Eqs. © and 
(0) we find 



\T rE {u(p- P >)m}\\ 2 ) = 



d\dE 



ds 



2(44 
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C(p), 



where C (p) denotes the generalized concurrence of p [23] . 
This shows that the average distance between the re- 
duced system states is proportional to the concurrence, 
i.e., to the amount of entanglement in the initial state. 
We note that Eqs. © and (TTTfl) are not intended to be 
used in actual experiments to determine 5(p), because 
the left-hand sides of these equations involve an average 
over all unitaries which is hard to realize experimentally 
The result © for the average increase of the Hilbert- 
Schmidt distance between open system states can be gen- 
eralized to other physically relevant random matrix en- 
sembles describing generic features of complex quantum 
systems. To this end, we write the unitary time evolution 
operator Ut — exp{— iHt} of the total system in terms 
of a unitary matrix W , formed by the eigenvectors of the 
Hamiltonian H, and a diagonal matrix D, containing the 



eigenvalues Ej of H, as Ut = W exp{— iDt}W'. We as- 
sume that the eigenvectors of H are random such that 
W becomes a random unitary matrix distributed accord- 
ing to the Haar measure. The characteristic behavior of 
a generic complex quantum system is determined by the 
distribution of its level spacing (see, e.g., [13]). In the fol- 
lowing we determine the average over W, working with a 
general level spacing distribution which includes, in par- 
ticular, the cases of regular and chaotic level statistics. 

Our aim is to determine the quantity on the left-hand 
side of Eq. © with U replaced by W exp{-iZM} . The 
averaging with respect to the Haar measure now requires 
the determination of an 8th moment of W. The calcula- 
tions can be performed efficiently by employing the tech- 
niques of Ref. [25[ which are based on advanced group 
theoretic methods and the Schur-Weyl duality. The re- 
sult may be written in the form 



|Tr E {We- iDt WHp - p')We iDt W^}\\ ) = cS(p), 

(11) 

where the prefactor c is independent of p, and depends 
only on the dimensions ds and dE, and on time t via the 



Fourier transform f(t) — 4^ . 



,-iE,-t 



of the level density. 



We conclude that, independently of the level spacing dis- 
tribution, the average Hilbert-Schmidt distance between 
the open system states is proportional to the measure 
for quantum correlations given by Eq. Hence, also 
for more realistic random matrix ensembles the reduced 
system distance increases on average if and only if there 
are initial quantum correlations. The features of specific 
matrix ensembles, describing e.g. regular or chaotic level 
statistics, can be investigated in further detail by deter- 
mining the corresponding averages of the function f(t). 
Recently, a similar approach has been pursued for the in- 
vestigation of non-Markovian complex open systems [2(| . 

In summary, we have developed a method to locally de- 
tect initial correlations in an open quantum system. The 
witness is successful on average if and only if the correla- 
tions are of nonclassical nature. This result nicely com- 
plements the recent results of Ref. [l6[ which show that a 
completely positive dynamical map describing the evolu- 
tion of the reduced system can be found if and only if the 
total initial state has zero quantum discord. The method 
developed here is based on non-selective measurements 
and state tomography of the reduced system at differ- 
ent times, and should thus be realizable with present-day 
experimental technologies [l8|, [lj| . Finally, we mention 
that the results of this work are not restricted to an open 
system setup, the method is applicable to any bipartite 
system. This includes typical communication protocols, 
where quantum states are shared by two parties. The 
correlations of the state can be detected locally by either 
one of the parties with the aid of the proposed procedure. 

Proof of the theorem. Let M = rrii \i) (i\ be the 
spectral decomposition of M, and an orthonormal 
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basis of the open system's Hilbert space Hs- The matrix 
elements A^i = ((pk\A\(pi) can be expressed as 

&id = 1 52mi(i\U i (\<Pi)(<Pk\®I)U\i), (12) 



which yields for the average Hilbcrt-Schmidt distance 



^l| a ) = EE m * m i< i l A A w At l (|i)0 , |) li)> (13) 



k,l i,j 

where we have introduced the map 



Kab{X)= \ dp{U)U^AUXU^BU, (14) 



and defined the operators A^i 
Lemma. We have 



\<Pl){<fik\ ® I- 



Aab(X) = a(TrX)I + bX 1 



(15) 



where 



d TrB A — TrA TiB , d Tr A TrB - TrB A 
a = 77-^ - , b = -rj-^ — . (16) 



d{d 2 - 1) 



d(d 2 - 1) 



Although the proof of this lemma can be carried out 
with the help of the general method of Ref. [25j ]. here 
we present an elementary proof. To this end, we use the 
unitary invariance of the Haar measure to show that the 
map Aab has the invariance property Aab{W XW*) — 
WAab(X)W^ for all unitary operators W. Under the 
Choi- Jamiolkowski isomorphism [27], [28| Aab is mapped 
to the operator p\ AB = {Aab ® I) where |fi) = 

\co) ® \u}) £ H ® H is a fixed maximally entan- 



1 ^ d 



s/d 



gled state. The above unitary invariance of Aab implies 
that pa ab is U ® {/""-invariant or isotropic 2t| 3pj . with 



the complex conjugation defined with respect to the ba- 
sis The general form of such operators is known 
to be a linear combination of the unit operator and the 
projection onto the state p\ AB = ol/d + b |fi)(fi|. 
The constants a and b can be obtained by determining 
the trace of pA AB and the matrix element (n\pA AB \Q) . 
Applying the inverse Choi- Jamiolkowski isomorphism to 
PA AB , one finds the general form (I15|) for the map Aab, 
with a and b given by Eq. (|16[) . which proves the lemma. 

We now insert the expression (TT51) into Eq. (fT3")l . Using 
the relations Tr{A^A fc j} = d E and TrA fei = TrA^ = 
SkidE one can see that the resulting expression reduces 
to Eq. ([8]), which concludes the proof of the theorem. 
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